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A hierarchy of commutative Poisson subalgebras for the Sklyanin bracket is proposed. Each 
of the subalgebras provides a complete set of integrals in involution with respect to the Sklyanin 
bracket. Using different representations of the bracket, we find some integrable models and a 
separation of variables for them. The models obtained are deformations of known integrable systems 
like the Goryachev-Chaplygin top, the Toda lattice and the Heisenberg model. 

1 Introduction. 

Let us consider a 2 x 2 matrix T(A) which depends polynomially on the parameter A 

T(A) (1-1) 

whose entries are polynomials of the form 

A(X) = X N + A N ^X N - 1 + ... + A , B (A) = B N ^ X N ^ + . . . + B , (1.2) 



C (A) = C^-i A"" 1 + . . . + C , D (A) = D N _ 2 X N ~ 2 + . . . + D . 

The algebra 21 at of all functions depending on the coefficients 

A ,...,A N _i, B ,...,B N _ 1 , Co,...,Cn-i, Dq,...,D n _ 2 - (1-3) 

of the matrix T(A) is a Poisson algebra with respect to the so called Sklyanin bracket ||. 
The dimension of generic simplectic leaves for this bracket is equal to 2N. 
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The main property of the Sklyanin bracket is that the coefficients of the trace of T(A) 
generate an iV-dimensional commutative subalgebra Ql° N in 21 jy The generators of 2ljy are 
linear polynomials. In this paper we construct different iV-dimensional commutative subal- 
gebras 2tjjf , M 6 N, in 21jv generated by polynomials of higher degrees. These subalgebras 
can be regarded as deformations of the standard trace subalgebra 2l^r • The generators of 
Sljlf define a set = {If 1 , • • • , Ijf} of N integrals in involution. The separation variables 
for all these sets S^f coincide with those for S% whereas the separated curves are different. 

Using known representations for the Sklyanin bracket and our families Sf/ of integrals 
in involution we can construct new integrable models. They are deformations of integrable 
models corresponding to the standard trace integrals. One of the examples is related to 
known representation 0, [L| of 2I2 on the Poisson algebra e(3). 

The most interesting models are related to the commutative subalgebras 2t]y possessing 
a quadratic polynomial in variables ( |1.3|) which can be regarded as a Hamiltonian. For 
example, the subalgebra Oil gives us the quadratic deformation of the Goryachev-Chaplygin 
top found in 0. In this paper we construct integrable deformations of arbitrary degree for 
the Goryachev-Chaplygin Hamiltonian. 

The well-known Goryachev-Chaplygin case (with the additional gyrostatic term) in rigid 
body dynamics is described by the following Hamiltonian 

H = J 2 + J 2 + 4 J| + 2 Cl x x + 2c 2 x 2 + c 3 J 3 , (1.4) 

where q are arbitrary constants. The Lie-Poisson brackets 



{Ji j Jj } £-ijk Jk 1 \Ji j %j } £-ijk %k 1 1 -Ej } 0; ^ j ji k 1)2,3, (1.5) 

where is the standard totally skew-symmetric tensor, defines the corresponding equa- 
tions of motion. These brackets possess two Casimir elements (x, x) and (x,J), where 
J = (Ji, J2, J3), x — (xi, X2, X3) and (x, y) stands for the scalar product in M 3 . 

On the fixed level (x, J) = of the second Casimir element the Hamiltonian (|1.4j ) 
commutes with an additional cubic integral of motion. This fact ensures the integrability of 
the Goryachev-Chaplygin case. 

The quadratic deformation of the Goryachev-Chaplygin Hamiltonian given by 

H = h = Jl + Jf + 4 J 3 2 + 2 Cl x x + 2c 2 x 2 + c 3 J 3 + 

(1.6) 

4(aiXi + a 2 x 2 ) J 3 - (a? + a 2 2 )x\ 

also has an additional cubic integral [0]. If C\ = c 2 = c 3 = then ( |1.6|) gives us a new 
partially integrable case (i.e. integrable on a special level (x, J) = of one of the integrals 
of motion) for the Kirchhoff problem of motion of a rigid body in the ideal fluid. A similar 
deformation of the Kowalewski top has been considered in || [jj. 

As an application of our general scheme we present below a separation of variables for 
this model. The canonical separated variables for the Goryachev-Chaplygin top @ are given 
by 



q 1>2 = J 3 ± J Jf + Jf + Jf , p h2 = —In (q lj2 (ixi - x 2 ) - (iJ 1 - J 2 )x 3 ) . (1.7) 
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If we substitute two pairs of variables (|1.7|) into the following separated equation 

A 3 + a A 2 - h\ + I = c fi — , A = q lj2 , fi = exp(2ip h2 ) , (1.8) 

and solve the couple of linear equations obtained with respect to I± = H and Jo we immedi- 
ately derive the Hamilton function 

H gch = J\ + Jf + 4J 3 2 + 2a J 3 - i{c + M^i + (c - b )x 2 • (1-9) 

and additional integrals of motion for ( |1.9| ). 

Substituting the same variables (p, q) into another separated equation 

A 3 + a A 2 - h\ + I = ( Cl A + c )/i - (M + ^ { ^ %) , (1.10) 

H 

we get an integrable system with the following Hamiltonian 

H = H gch + (ci(Ui - J 2 ) + bi(Ui + J 2 )^jx 3 - 2(c 1 (ix 1 - x 2 ) + b^ixi + x 2 j) J 3 . (1.11) 
After the canonical transformation 

x^x, J^J + Ux, U= I -c_ I, c± = bl ^ C \ (1.12) 



the latter Hamiltonian becomes 









— ic + 








— c_ 


ic + 


c_ 





)Xi 


- (ci 


-h) 



-i(b + c - a (ci + bi))x! + (c - b - a (ci - &i))x 2 , 



(1.13) 



which coincides with (|T! 

In Section 4 the simplest deformations for the Toda lattice and the Heisenberg model 
are given. We do not know whether physical applications of such deformed models exist. 

Acknowledgments. The authors are grateful to the Newton Institute (Univ. of Cam- 
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and 99-01-00698, INTAS grants 99-1782 and 99-01459, and EPSRC grant GR K99015. 



2 Properties of the Sklyanin bracket 

The explicit form of the Sklyanin brackets for the coefficients of the matrix (|1 . 1| ) can be 
derived from the following operator definition 

{T(A), T{y)) = [r(A - a*), r(A)r(A0 ] , (2.1) 
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1 2 

where we use the standard notations T(A) = T(A) ® Id , T(/i) = Id® T(/i) and 



r(A - /j) 



7/ 



A — /i 



n 



where 



n 



/ 1 o o o \ 

10 

10 

\ 1 J 



t] eC 



(2.2) 



Example 1. In the simplest case N = 2 relation (2.1) is equivalent to 

{A 1 ,A Q } = 0, {A 1 ,B 1 } = r ] B 1 , {A 1 ,C 1 } = -rjC 1 , {A u B } = rjB , 
{A 1 ,C } = -T}Co, {A 1 ,D o } = 0, {A ,B 1 } = r ] B , {A , C x } = -r}C , 
{A , B } = r l (B A 1 - A B 1 ), {A , C } = ?,(AA - C A 1 ), (2.3) 
{A , D } = ^BoC, - fliCb), {B u d} = {5!, B } = 0, 

{Si, Co} = -rjDo, {B u D } = 0, {Bo, C x } = - V D , {B , C } = -r ] D A 1 , 
{Bo, D } = -vB.Do, {C x , Co} = {C u D } = 0, {C , D } = rjC^. 



It was proven in (0) that the coefficients of the determinant 
d{\) = detT(A) = A(X)D(X) - B(X)C(X) 



(2.4) 



belong to the centre of 21 or, in other words, they are Casimir elements for bracket ( |2.1| ). The 
number of the Casimir functions is 2N—1 and therefore we have a AN— 1 dimensional Poisson 
manifold with degenerate Poisson structure ( j2.1| ) and 2iV-dimensional generic symplectic 
leaves. 

To bring the Poisson bracket ( |2.1| ) into canonical form, a new set of variables 

d , ■ ■ • , d 2 N-2, qi,---q N , pt,---,p N . (2.5) 
was proposed in The variables do, ■ ■ ■ , d 2 N-2 are the coefficients of fl2_. 

d(X) = d 2N ^ 2 X 2N - 2 + ... + d 
and the variables qi are zeros of the polynomial A(X): 

N 

A(X) = H(X- qj ); 

i=i 

Pj = V In-B(^) . 



the variables Pi are defined by 
It follows from (O) that 



{A{\),A(ii)} = {B{\),B(ri} = 



and 



{A(X),B(»)} 



X — jJL 



A{X)B^)-A^)B{X)). 



These relations imply {qj,qk} = {Pj,Pk} = and {pj,qk} = $jk j respectively. As usual if 
we fix values of the Casimir elements dj, then the canonically conjugated variables qi,Pi are 
simplectic coordinates on the corresponding simplectic leaf. 

To express the variables ( O ) in terms of (f2.5|), one can use the formulae 



3 A construction of commutative subalgebras 

Let us introduce the matrix 

T(X) = K(X)T(X), (3.1) 

where T(X) is given by fll.lj) , 

*(A) = ( c ) (A) ' (3 ' 2) 
whose entries are polynomials of the form 

A (A) = A M X M + Am-iX* 1 - 1 + ■ ■ ■ + Ao , 

B (A) = b M X M + bM^X" 1 - 1 + ■ ■ ■ + b , C (A) = c M X M + c M -iX M - 1 + ■ ■ • + c , 

and bk, Ck G C are arbitrary parameters. We require that the trace of the new matrix T has 
the form 

N+M N-l 

trace f (A) = ^ a^X 1 + ^ J fc A fc , (3.3) 

i=N k=0 

where a, are arbitrary fixed constant parameters. It easy to see that both unknown coeffi- 
cients Ai of the polynomial A and unknown functions Ik in ( |3.3| ) are uniquely defined from 
( p.3|) . Moreover, they are polynomials in variables (|1.3| ) such that all the polynomials Ik 
have the same degree M + 1 and the degree of equals M — i. 

Example 1 (continuation). If N = 2 and M = 1 the functions .Ai and Jj are given 

by 

Ax = ai, Ao = a — d\A\ — b\C\ — ciBi , 

h = ax(A - A\) + 6 a (C - Axd) + ci(B - MBx) + aoA + 6 Ci + c 5i , (3.4) 

= (o - a-iM - b x Ci - cxBi)Ao + c B + b C . (3.5) 

Theorem 1. Polynomials Jj defined by Q3.3|) commute with each other with respect to 
the bracket (12.11). 
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Proof. The explicit form of condition ( |3.3|) is given by 

N+M 



N-l 



A(X)A(X) + B(X)C(X) + C(X)B(X) = a i-M^ + Yl /fcA "'- 



k=0 



Let us substitute A = qj into this identity. It follows from the definition d(X) = AD — BC 
that 

C{qj) = = -d(Qj) ex P(-Wj) • 

Taking this formula into account we get TV linear equations 

N-l N+M 

C(qj) exp(?/p i ) - B(qj) d(qj) exp(-ljpj) = ^ + Y a k-N Qj , j = l,...,N. (3.6) 



8=0 



k=N 



for TV unknowns 7j. Following ||, we can easily show that the functions 7j are in involution 
with respect to Q2.1|) . 

Let us rewrite this system in the matrix form 



where 



and 



<f>(p,q) = S(q)I(p,q) + U(q), (3.7) 

C(qj) exp(r] Pj ) - B{qj) d{ qj ) exp(-r] Pj ) , Uj = ^ a k Qj 

k 

/i, gi, ••• 



\ 1, Qn, ■■■ Qn 1 J 
The solution of ( |3.7| ) is given by 

J(p, q) = W(g) ( 0(p, g) - 17(g) ) , W(g) = S" 1 ^). (3.8) 

The matrix S belongs to the class of so-called Stackel matrices || for which the entries 
Sij of the z-th row are functions depending on g, only. Differentiating the identity 

W mi Su = 5 m i , 

i 

with respect to qj and using the main property of the Stackel matrices, we have 

Substituting k for m and eliminating the right hand side dSji/dqj from the two obtained 
equations, we find 
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Since the determinant of the matrix S does not vanish 



„, dW im dW ik 



for all j, k, i, m. 

Calculating the Poisson brackets between I k and I m , we obtain 

{'..'»> - £ - wj^) f*(* - UM) • (3-10) 

It follows from (|3.9| ) that I m } — . The proof is complete. 

Collorary 1. Variables Pi,qi are the separated variables for integrals of motion Ij 
satisfying the separated equations (|3.6| ) -(|3.7| ) whereas variables 

\i = q h fii = C(qi) exp (r)Pi) 

satisfy the characteristic equation Det(T(X) — jj) — 0. This means that these variables lie 
on the corresponding equivalent algebraic curves. 

Remark 1. Theorem 1 for M = is a well-known fact. In the case N = 2, M = 1 the 
matrix T with 

X + Ao b 



K ^ = V c ) ' A = a Q -A 1 , (3.11) 

has been proposed in in a non-factorized form. In the paper Q this matrix was factorized 
and generalized to the case of arbitrary N. Notice, that our matrix K for M = 1 has a more 
general form 

K(X) = ( * M + bo Y A Q = a -A 1 - hd - Cl B u 



than ( |3.11| ). Parameters b\ and c\ are absolutely essential in constructing new quadratic 
integrable Hamiltonians. 

Remark 2. To generalize the condition ( |3.3D we can assume that 

trace T(A) = ^a k X k + ^1^, a k G C . (3.12) 

Here g = {ii, . . . , i]y} is an arbitrary subset in the set of numbers {0,1, . . . , N + M} and g 
is the corresponding complement such that g U g = {0, 1, . . . , N + M}. But in this case ij 
are rational functions of variables ( |1.3| ). 

Example 1 (continuation). If N = 2 and M = 1 and 

trace T( A) = a! A 3 + I\X 2 + a A + I , 

these functions are 

Ai = ax, Aq = ^- (a - aiA - cqBx - boCx - biC - ctBo^j , (3.13) 



/i = J- i ( a °- a Mo ~ A\) - h(C - dAx) - ci(S - #i^i) - co^! - bod] 
Io = J-^Mao - cxBq - ai A - hC ) + &o(C Ai - 4A) + c (5 ^i - A £i)) . 



Remark 3. Theorem 1 can be proven in the same way for matrices (|1 . 1| ) of slightly 
different structure. For example, we can assume that the entry D(X) of the matrix T(A) has 
the form 

D(X) = X N + D^X"- 1 + D N _ 2 X N ~ 2 + • • • + A> . 
In this case we have one additional variable -Djv-i and one additional Casimir function 

d-2N-l = An-1 + Dn-1- 



Polynomial deformations of known integrable mod- 
els 



If we identify the Sklyanin bracket ( |2.1| ) with a fixed Poisson bracket on some phase space 
A4, then ([2.1|) can be regarded as an equation for matrix T(A) of the form (|1.1|) whose 
entries ( |1.2| ) are polynomials in A with coefficients being functions on M.. Such a matrix 
T(A) defines a representation of (|2.1|) on M.. 

Using known representations and Theorem 1, we can produce hierarchies of deformations 
for several integrable models such as the Goryachev-Chaplygin top, the Toda lattice, and 
the Heisenberg magnetic. According to Collorary 1 integrals of motion for all members of 
such a hierarchy are separable in the same canonical variables Pi,qi- 



4.1 Deformed Goryachev-Chaplygin top 

Let us consider the Lie algebra e*(3) with a natural Lie- Poisson bracket (|1.5| ). The phase 
space of the Goryachev -Chaplygin top M. is a union of special non-generic coadjoint orbits 
(symplectic leaves) of E(3) in e*(3) specified by the fixed value I 2 = of the second Casimir 
operator. 

It was observed in |I[ that the functions 

Ax = -2J 3 , A = -J\ - Jf J; B L = x 2 + ixx, B = -(J 2 + iJi)x 3 , 



2 ' 
X\ 

Cx = —x 2 + ixx, C = (J2 - iJi)xs, D = x\ 



satisfy ( p.3| ) if the Poisson bracket is given by ( |1.5| ), I 2 = and 77 = —2i. Formulae 
and (|3.5|) give us two integrals Ji and J such that {ii, J } = 0. It easy to verify that (up to 
constant factors) 



Hx = h = ax [Jx + Ji + ^1 + ~2j + 2ao^3 + b (x 2 - ix x ) - c Q (x 2 + ixx) 
+bx(2J 3 x 2 - J 2 x 3 - 2iJ 3 xx + iJxx 3 ) - cx(2J 3 x 2 - J 2 x 3 + 2iJ 3 x x - iJxx 3 ) 



(4.1) 



which is equivalent to ( |1.6|) . The explicit form of the integral Io is given by 

Io = 2ax(Jf + Jf) J 3 + ao{Jl + J|) + b (-J 2 + Ux)x 3 + c (J 2 + Ui)x 3 
+bx{Jx + J 2 )(x 2 - ixx) - Cx{Jl + J 2 ){x 2 + ixx) 
2axJ 3 + a + bx(x 2 - ixx) - Cx(x 2 + ixx) 

I £ o * 

x 3 



(4.2) 



As examples we also present explicitly one polynomial deformation of higher degree and 
one rational deformation. Thus if M = 2 the deformation reads as follows: 

H 2 = H 1 +Aa 2 Jz({Jl + Jl + 2Jl)+exf^ 

+ b 2 (2iJ 3 x 3 {Jx + U 2 ) + (Ji + J 2 + 4J| + ex 3 2 ){x 2 - ixif) 

- c 2 (^2U 3 x 3 (J 1 - iJ 2 ) + (J 2 + J 2 + 4J 3 2 + exf)(x 2 + ix x f), 
where H\ is defined by ( |4.1|) . In the case M = 1 the simplest rational deformation 

g = Hi- «o(2J 3 - 1) 



2J 3 

corresponds to (|3.13|) . 

4.2 Deformed Toda lattices 

The Sklyanin bracket may be identified || with the standard bracket in M. = H 2N with the 
help of the following ansatz for the matrix T(A): 

T(A) = L N Ljy-i ■■■L 1 , Li(X) = ^ X _ e - q l 6 q J Vi+N = Qi- ( 4 -3) 

Here Pi, qi are canonical variables in M? N . 

The matrix T(A) obeys the Sklyanin bracket ( |2.1| ) with r\ — 1 and describes the periodical 
Toda lattice 0. The trace of the matrix T(A) produces the integrals of motion. The simplest 
of them are 

N N N 

P = -Y^Pi, H = -J2 PiPj -J2e qi+1 - qi ■ (4.4) 

i=l i>j i=l 

The matrix ( f4.3|) has the necessary polynomial structure ( |1.2|) and we may apply Theorem 
1 in order to produce deformations of the periodic Toda lattice, polynomial in momenta. 
For instance, if M = 1 then the integral of the lowest degree in trace T(A) is the following 
quadratic integral 

N N N 

Jjv-i = -aoJ^Pi ~ oi(X)P*Pi + J2 eqi+1 ~ qi ) + e91 ( c o + c iPi) - e ' QN ( b o + b iPN) • (4.5) 

i=l i>j i=l 

4.3 Deformed spin chain 

Let the phase space M. be a direct sum of the algebras si* (2) with the following Lie-Poisson 
bracket 

{4, 4} = ±4, = 24 (4.6) 
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in the each algebra. The Sklyanin bracket with rj — 1 is identified with ( |4.6| ) by 



T(A) — L N L N _i ■ ■ ■ Li , Li{\) 



A + 4 
si A 



It is easy to check that such a product has the polynomial structure described in Remark 
3 under the restriction that the Casimir fumction d>2 n -i = A/v-i + Dn+i is equal to zero. 
Adapting our construction to this case we can construct polynomial deformations for the spin 
chain. In the simplest case M = 1 we obtain the following integrable quadratic Hamiltonian 



TV 

i \ - i — H — a S 3 + b S + + c S^ + ai (^2(s l + s j _ + s l 3 s J 3 ) - SfJ 

j>i 

N . N 



j>i i=l 
N ^ N 



j<i i=l 



where 



S3 = An-i = ^ S3, S- = -Btv-i = s % _, S + = Cn-i = ^ s\. 



If bi = cx = then H coincides with the general quadratic integral for the standard spin 
chain. 

It is well-known that the Hamiltonians for the Toda lattice and the spin chain derived 
from the factorised form of operator T(A) have to be transformed to bring them to a local 
form suitable for applications. We do not know such transformations for the deformed 
models. 
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